The Phase Diagram of the Shastry-Sutherland Antiferromagnet 
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The Shastry-Sutherland model, which provides a representation of the magnetic properties of 
SrCu2 (B03)2 , has been studied by series expansion methods at T = 0. Our results support the 
existence of an intermediate phase between the Neel long-range ordered phase and the short-range 
dimer phase. They provide strong evidence against the existence of helical order in the intermediate 
phase, and and somewhat less strong evidence against a plaquette-singlet phase. The nature of the 
intermediate phase thus remains elusive. It appears to be gapless, or nearly so. 

PACS numbers: 75.10.Jm., 75.40.Gb 



I. INTRODUCTION 

The recent discovery of the two-dimensional spin gap system SrCu2(B03)2 [Pj, and its representation as an antifer- 
romagnetic S = \ Heisenberg spin model |2| equivalent to a model previously introduced by Shastry and Sutherland 
[pf , has led to many recent studies of this model. 

The Shastry-Sutherland model is a nearest-neighbor square lattice antiferromagnet, with additional diagonal inter- 
actions, in a staggered pattern, on alternate squares. The Hamiltonian is written as 



H 



diag 



(1) 



axial 



This convention is adopted because in SrCu2 (603)2 the diagonal bonds are, in fact, the shortest and the exchange 
constants are estimated to be J ~ 100K , J' ~ 68K. The model is illustrated in Figure [[} 
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FIG. 1. The Shastry-Sutherland lattice. The exchange J' acts between sites separated by horizontal and vertical links, while 
the exchange J acts across the diagonal dashed links. The spin orientation at each site for helical order is also given (near the 
sites), where q = arccos(— J'/J). 



It is clear that the system will exhibit Neel order for J/ J' small, and will form a gapped spin-singlet state for J/J' 



large, with S = dimers on the diagonal bonds. Shastry and Sutherland showed that this dimer state is an exact 
eigenstate of the Hamiltonian for any J, J' and that it is rigorously the ground state for J/J'>2 |. 

The main controversial and challenging question is whether the system has an intermediate phase, and what kind 
of phase it is, if it exists. In the classical limit, when S — > oo, the ground state is Neel ordered if J/ J' < 1 and 
is helically ordered (see Fig. ft]) otherwise, where the twist between one spin and its nearest neighbor is given by 
q = arccos(— J'/J). For the S = \ system, Albrecht and Mila Q used a Schwinger boson mean-field theory (SBMFT) 
and found an intermediate phase with helical long-range order (LRO) for 1.1 < J/ J' < 1.65, but with q differing from 
its classical value. In our previous series study 0, we were unaware of this work, and did not consider the possibility 
of a helical LRO intermediate phase. We computed series for the Neel-ordered phase and for the dimer phase at 
T = 0, and located a transition point at J/ J' — 1.45(1) which appeared to be first-order. 

Subsequently Koga and Kawakami |6J calculated series expansions about disconnected plaquettes, and claimed to 
have identified two transition points: a second-order transition from Neel order to a plaquette singlet phase at J/ J' ~ 
1.16, followed by a first-order transition to the dimer phase at J/ J' ~ 1.48. We find this interpretation unconvincing, 
partly because of uncertainty in the analysis of rather short series, and partly because such an intermediate phase 
is totally different from the helical-ordered phase presented in Ref. ||]. Another recent development is the field- 
theoretical study of a generalized model with Sp(2N) symmetry by Chung, Marston and Sachdev [Q. For S = |j 
they suggest a Neel to helical LRO phase transition at J/ J' ~ 1.02, and a helical LRO to dimer short-range ordered 
(SRO) phase transition at J/ J' ~ 2.7. This is similar to the scenario in Ref. 0], but with the helical phase extending 
over a much larger range. The quantitative accuracy of this approach for S = ^ is problematic. Both our previous 
work p| and the plaquette series of Koga and Kawakami jq| provide strong evidence for the dimer phase setting in at 
around J/ J' ~ 1.45 — 1.48. It is also interesting to note that for 1/S > 5 (S is a continuous variable in this theory) a 
phase with plaquette SRO is predicted H] . If finite N fluctuations change the phase diagram significantly then it is 
conceivable, according to this theory, that there are four stable phases. 

Knetter et al. [|3| subsequently performed a new series analysis of the dimer phase. They find that the lowest 
single-triplet excitation vanishes at J'/J — 0.69 (or J / J' = 1.45), very close to the transition point mentioned above; 
but that an S = 1 two-triplet excitation energy vanishes even earlier, at J'/J = 0.63 (or J/ J' — 1.59). This would 
seem to indicate important binding between the triplets, possibly leading to a condensate of triplets at or even below 
J'/J — 0.63. Even more recently, Totsuka, Miyahara, and Ueda || have discussed two-triplet binding effects in this 
model using both perturbation theory and exact diagonalization, and have shown that binding occurs even in the 
quintet channel of two triplets. 

All the above studies show that the quantum system is Neel ordered up to and beyond J/ J' = 1, i.e. beyond the 
regime for the classical system: this is also consistent with the general phenomenon of quantum order by disorder 
JUM. The transition to dimer order is located at J/ J' ~ 1.45 — 1.58, so if there is an intermediate phase, it must exist 
within the region 1 < J/ J' < 1.5. Our aim, in the present paper, is to explore these issues through new extended 
series expansions. We have obtained extended series expansions about isolated plaquettes, both with and without 
diagonal bonds, following Koga and Kawakami M. Secondly, we have extended our previous series about the Neel 
ordered state and the dimer singlet state H . Finally, we present new series expansions about states with helical LRO 



and columnar dimer order. The next Section presents technical details of the calculations, with analysis of series and 
discussion of results given in Section III. In the last section we present our conclusions. 

We conclude that an intermediate phase between the Neel phase and the dimer phase very likely exists, in the 
region 1.2 < J/ J' < 1.5. The nature of this phase is much less certain, however. It does not appear to be plaquette 
ordered or helically ordered. It could perhaps posses weak columnar dimer order; or it could posses some other type 
of order, not considered here. In any case, this intermediate phase appears to be gapless, or nearly so. 

II. SERIES EXPANSIONS 

We use the linked-cluster expansion method to derive perturbation expansions for various choices of the unperturbed 
reference Hamiltonian. The technical details are discussed in many articles, and we refer the reader who is unfamiliar 
with these to a recent review fllfl . The following is a brief summary, 

A. Dimer Expansions 

The first term in the Hamiltonian (|lj), which represents disconnected dimers, is taken as the unperturbed Hamilto- 
nian. The unperturbed ground state is then a product state of S = singlets on these dimers. As mentioned above, 
this state remains an exact eigenstate of the system for all J, J', but is not the true ground state for J' /J > (J'/J) c , 
The perturbation, being the second term in (|lj), mixes this state with states in which triplet excitations occur on the 
dimers. 

Dimer expansions can be developed for both ground state properties and for excitations. Here because the dimer 
state is the exact ground state, we focus on the triplet excitation spectrum A(k). We have computed expansions to 
order (J'/ J) 21 , extending our previous calculation ||] by 6 terms. This calculation involves 185332 linked clusters up 
to 11 sites. The resulting series coefficients are available on request. The minimum gap is at k = 0, where the series 
is: 



A(k = 0)/J= 1 - x 



2 



X 



5x 5 7x 6 2051 x 7 39091 x 8 



2 8 32 384 4608 55296 

268849 x 9 964411a; 10 



2.180975296 a; 13 



663552 6635520 58982400 191102976000 

-1.938901500 x 14 - 0.071485040 a; 15 + 0.820962395 x 16 - 2.236236507 x 17 

-7.855506946 x 18 - 9.452785582 x 19 - 2.833667089 x 20 + 4.387303805 x 21 + (9(x 22 ) (2) 

where x = J'/ J. A standard Dlog Pade approximant analysis [1J] gives the gap vanishing at J' /J = 0.6965(15), or 
J I J' — 1.436(3) with the estimated critical index v ~ 0.45(2). Different Dlog Pade approximants show remarkable 
consistency with this exponent value, suggesting that this transition may belong to a new universality class. These 
results are in agreement with Muller-Hartmann et al. Jl3[ , who used shorter series. But we recall that Knetter et al. 
H found that a two-triplet 5 = 1 bound-state engery vanishes even before J' /J — 0.70, at around 0.63, which would 
entirely alter the position and the nature of this transition. We have not computed series for the two-particle states 
in this model. It would be very interesting to explore this transition further, using series or other methods. 



We also computed a new series, to order (J'/ J) 17 , for the susceptibility \ corresponding to this momentum k = 0. 
The resulting series is 

x 2 x 3 41 x 4 281 x 5 6461a; 6 65953 x 7 186863 x 8 

x(k = o)= i + -j + T + — 6 



1152 



13824 



165888 



311040 



374542669 a; 9 24573240371a; 10 2138459511091a; 11 „„„„„„„,„ 12 

+ — + 1.308076648 a; 12 



597196800 



28665446400 



2149908480000 



-1.604414343 a; 13 + 2.086275141 a; 14 + 2.638378225 x 15 + 3.431256777 x 16 + 4.419063264 a; 17 + 0(a 



(3) 



A standard Dlog Pade analysis of this series indicates a preferred critical point around J' /J — 0.71. Biasing x c at 
0.6965, the critical exponent is found to be very small, around 0.03, which is completely different from that of the 
classical 3D Heisenberg model. It might even be compatible with a logarithmic divergence. 

B. Plaquette Expansions 



(a) 



— 


\ ' 

\ 
\ 

1 ^ 


— 


t 1 

\ 
\ 

i ^ 


— 


T— 1 

— — 

/ 
/ 

/ 


I ' 

\ 

\ 
\ 





\ ' 


• 



(b) 



• ' 1 

\ 


i — i 


r 

% 

N 


• ( 


II i 






/ 
/ 












/ 
/ 
/ 






: k 




^ 


• i 



FIG. 2. (a) The first plaquette structure (PE1): The bold solid, the thin solid and the dashed lines represent the coupling 
constants J', A J' and A J, respectively, (b) The second plaquette structure (PE2): The bold solid, bold dashed, the thin solid, 
and the thin dashed lines represent the coupling constants J', J, XJ' and XJ, respectively. 



Instead of perturbing about isolated dimers, any set of disconnected units can be used. Using a plaquette basis 
allows the investigation of plaquette type order. Because each plaquette has 16 rather than 4 states it is not possible 
to derive series of the same length as with dimer expansions. 

Following Koga and Kawakami [p|, we have computed two kinds of plaquette expansions (PE1 and PE2), which 
respectively take the set of plaquettes without the diagonal J bonds and the set with diagonal J bonds as the 
unperturbed Hamiltonian. To make the expansion possible it is necessary to introduce an expansion parameter A < 1 
which modifies the interactions not included in Hq- This is illustrated in Figure 2. The series are then computed in 
powers of A. The analysis evaluates these at A = 1, corresponding to the original Hamiltonian. 

For PE1 we have computed series to order A 8 , A 7 , A 6 for the ground state energy, triplet excitation energies A(k) 
and staggered susceptibility xaf respectively. This is one additional term for Eq and two additional terms for A(k) 
and xaf over Ref. 0. The calculation is computationally demanding, even with our efficient program. For example, 
the computation of A(k) to order A 7 took about 10 days and required 1.8GB memory on an SGI Origin 2400 system 
with a 400MHz R12000 CPU. The effort for each additional term requires factors of approximate 50 and 15 increase 



in CPU time and memory, respectively. For the second expansion (PE2) the series have been computed to order A 7 
for the ground state energy and the singlet excitation spectrum, and to order A 6 for the triplet excitation A(k) and 
staggered susceptibility xaf- In Tables | and || we present the series for the ground state energy, the triplet gap A at 
k = and staggered susceptibility xaf for J/ J' = 1.25 and 1.4. Other series are available on request. The analysis 
is left for the next Section. 

C. Columnar Dimer Expansions 







Ik I 

\ 
\ 


* It II— -HI 

\ 
\ 

1 1- v - 




/ 
/ 




/ 
/ 

C i 


/ 
/ 






\ 

\ 


\ 






/ 
/ 




/ 
/ 
/ 1 


/ 
/ 






\ 
\ 

» — > 


\ 
1 II \ i> 



FIG. 3. Columnar dimerization pattern. The bold solid, the thin solid and the dashed lines represent the coupling constants 
J', A J' and A J, respectively. 

A plausible candidate for the intermediate phase is the columnar dimer phase, as shown in Fig. 3. For both the 
spin-i J i — J 2 Heisenberg model [|14| and the spin-i Heisenberg model on an anisotropic triangular lattice |lq| , there 
appears to be a columnar dimer ordered phase adjacent to the Neel ordered phase. To make the expansion possible 
it is necessary to introduce an expansion parameter A < 1 which modifies the interactions not included in Hq. This 
is illustrated in Figure ||. The series are then computed in powers of A. The analysis evaluates these at A = 1, 
corresponding to the original Hamiltonian. 

We have computed series to order A 10 , A 9 , A 9 for the ground state energy, the antiferromagnetic susceptibility Xaf, 
and the triplet excitation energies A(k) respectively. This calculation involves 290215 linked clusters up to 10 sites. 
In Table | we present the series for the ground state energy xaf f° r J/ J' = 1-25 and 1.4. Other series are available 
on request. The analysis is left for the next Section. 

D. Ising Expansions for Neel order 



For small J / J' , the model will have a Neel ordered phase with long range antiferromagnetic order (in the z direction). 
To construct an Ising expansion, we introduce an exchange anisotropy parameter A, and write the Hamiltonian as 



H = H + \V 



(4) 



where 

H = J 22 S i S j + J' z2 S t S J ' 

diag axial 

V = JJ2(S?S? + S»S») + J'J2 (S*S* + SfS*) , (5) 

diag axial 

The unperturbed ground state is the classical Neel state with energy Eq/N = —J'/2 + J/8. The perturbation V flips 
pairs of neighbouring spins. Series are developed in powers of A and evaluated at A = 1, which recovers the original 
Hamiltonian. We have obtained series for the ground state energy per site, Eq/N, and the staggered magnetization 
(order parameter) M to order A 12 , extending our previous calculation by 3 terms. We have also computed a new 
series, to order A 12 for the perpendicular susceptibility x±- The resulting series for the ground state energy and the 



staggered magnetization (order parameter) M for J / J' — 1.25 and 1.4 are listed in Tables || and [11; the series for 
other values of J / J' are available on request. For details of the analysis we refer to our previous paper ||. Results of 
the analysis are presented in the next section. 

E. Ising Expansions for Helical Order 

As discussed by Albrecht and Mila ||, the classical system has planar helical order for J /J' > 1. Starting from a 
reference spin in the z direction, each neighboring spin is rotated by an angle q, as shown in Figure 1. The twist is 
determined by minimization of the energy and yields 

fir, J < J'; 

9={ (6) 

[ a,rccos(-J'/J), J > J'. 

To develop an Ising expansion about such a helically ordered state we transform the Hamiltonian by rotating the 
spin axes at each site. The transformed Hamiltonian is 

H = H 1 + \(H 2 +H 3 ), (7) 

where 

Hi = Jcos(2q) J2 SISn + J' cos(q) Y, S i S j . 

diag axial 

H 2 = JJ2 [SfSl + co S (2q)S?SZ + sin(2q)(5f S* - S*S* n )} , (8) 

diag 



H 3 = J' J2 [Si Si + oos(q)SfSJ + sm(q)(SfSJ - S?S*)] . 

axial 

and where exchange anisotropy is introduced through the perturbation parameter A. 

We have computed series for the ground state energy and the order parameter to order A 12 , for various choices 
of q and the ratio J/J'. The resulting series for the ground state energy and the staggered magnetization (order 



parameter) M for J/ J' = 1.25 and 1.4 are listed in Tables [| and [II; other series can be supplied on request. We are 
not aware of previous series of this type for this model. 



III. RESULTS AND DISCUSSION 



In this Section we present a variety of results obtained from analysis of the various series expansions. The analysis 
has been carried out using integrated first-order inhomogeneous differential approximants and Pade approximants |l2J] 
to extrapolate each series to the physical value A = 1. 
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FIG. 4. The ground state energy versus J/ J' obtained from the Neel ordered Ising expansion (N), columnar dimer expansion 
(C), plaquette expansion (P), and helical order Ising expansion (H). The energy of the exact dimer state is also presented (D). 



A. Ground State Energy 



The nature of the ground state for any particular J/ J' can be determined by comparing the energies of various 
candidate states. Our analysis shows that, among all the series expansions we have calculated, the ground state 
energy from the Ising expansion about Neel order has the lowest energy for coupling J/ J' < 1.45, while the ground 
state energies obtained from the columnar dimer expansion and from the Neel order Ising expansion remain nearly 
equal through J/ J' > 1.2. For example, we estimate the ground state energy at J/ J' = 1.25 as 

-0.5628(15), Neel order; 

-0.5586(10), helical order; 

E /NJ' = { -0.5599(3), PE1; (9) 

-0.5584(8), PE2; 

-0.563(1), columnar dimer. 

and at J/ J' = 1.40 as 



(10) 



-0.561(3), Neel order; 

-0.554(1), helical order; 

E /NJ' = { -0.5548(5), PE1; 

-0.5525(10), PE2; 

—0.557(2), columnar dimer. 

where we have used the classical-g value in Eq. (ra) in the calculation of helical order. In Table IV we show details of 



the analysis from the latter case J/ J' = 1.4 for both the Ising expansion about Neel order and the plaquette expansion 
(PE1), from which we estimate the values in (UQ). The ground state energies obtained from various expansions versus 
J I J' are given in Fig. 0. Since the energy from the second plaquette expansion is slightly higher than that from 
the first plaquette expansion, we do not show the results of the second plaquette expansion. Our results are thus in 
disagreement with ref . H , which claimed that the plaquette phase has the lowest energy. Of course these energies are 
all very close to each other and the estimated errors are subjective confidence limits. However inspection of the data 



in Tabic IV shows that the Ising expansion about Neel order consistently gives lower energy estimates. 



Our results also argue against helical order being a stable ground state. In Figure H we show estimates of the ground 
state energy from Ising expansions about helical order, as a function of angle q, for the two coupling ratios J/ J' = 1.25, 
1.40. In each case the curve shows no indication of a minimum at some q c , but rather decreases monotonically towards 
q = ir, corresponding to Neel order. 
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FIG. 5. The ground state energy versus q for helical order. 



Over the region 1.2 < J/ J' < 1.45, the ground state energies obtained from the Neel ordered Ising expansion 
and from the columnar dimer expansion have some overlap after considering the error bars, so they are both good 
candidates for the ground state of the Shastry-Sutherland model. 



B. Staggered Magnetization and Perpendicular Susceptibility 

The staggered magnetization and perpendicular susceptibility will be non-zero in a phase with long range antifer- 
romagnetic order, and are expected to vanish at a transition point to a magnetically disordered or spin-liquid phase. 
Effective Lagrangian theory predicts a relationship p s = v 2 x± where p s is the spin stiffness and v the spin- wave 
velocity. In Figures o we show estimates of M and x± versus J/J', from the Ising expansion about Neel order. It 
can be seen that both quantities behave similarly, decreasing from their values at J = and vanishing at around 
J J J' = \.2 ± 0.1. The error bars are large and the vanishing point can not be obtained to high precision. However 
our new extended series shows a much sharper drop-off in M than given in our previous work B . Thus it is possible 
that the true transition point is at or below 1.20. 

For the Ising expansion about helical order the magnetization M at A = 1 (for the classical q value) is zero, within 
error bars, for all J / J' . This strengthens our conclusion, from the ground state energy results, that a helical phase is 
not present. 





FIG. 6. The staggered magnetization M and the perpendicular susceptibility xi versus J/J'. The solid points with error 
bars are the estimates from the Ising expansion about Neel order. 



C. Energy gap and susceptibility in the plaquette expansion 



The main argument of Koga and Kawakami H for the existence of a plaquette intermediate phase in the system 
is based on analysis of the triplet gap at k = and the staggered susceptibility series. In their series analysis, they 
assumed that the minimum triplet gap is at k = 0, and the transition lies in the same universality class as the classical 
3D Heisenberg model (i.e. critical exponents 7 ~ 1.4, v ~ 0.71). Using standard Dlog Pade approximants, they found 
an apparent critical singularity at A > 1 for J/ J' > 1.16. This would imply a non-zero spin gap for J/ J' > 1.16 
for the original Hamiltonian (hi). The above assumption is valid for the transition to Neel order, but if there was 



an intermediate phase, the transition from the plaquette phase to this intermediate phase might not lie in the same 
universality class as the classical 3D Heisenberg model, and the minimum triplet gap might not be at k = 0. In the 
analysis with our longer series, we find indeed that the minimum triplet gap is not at k = for J j J' > 1.25, and also 
that the transition does not lie in the classical 3D Heisenberg universality class. 
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FIG. 7. The triplet dispersion for various A and J / J' = 1.40 along the line connecting k = (0,0) and k = (tv/4,3tt/1&) 
obtained from the first plaquette expansion (PE1). 

Firstly, let us discuss the first plaquette expansion (PE1). Here for J j J' < 1.25, or for larger values of J j J' with 
small values of A, we find the minimum triplet gap is at k = 0; but for J / J' > 1.25 and A ~ 1, it is no longer at k = 0. 
For example for J j J' = 1.40, the dispersion for various A along the line connecting k = (0,0) and k = (w/4, 37r/16) 
is shown in Fig. [j], where we can see that for A < 0.75, the minimum gap is located at k = 0, while for larger A, the 
minimum gap is located at about k = (7r/4, 37r/16), and the dispersion near the minimum is quite flat. The series for 
the triplet gap at k = (7r/4, 37r/16) is 

A(tt/4, 3tt/16)/J'= 1 - 0.3077152787A - 0.2293662299A 2 - 0.0346064996A 3 

-0.0741227651A 4 - 0.0176489659A 5 - 0.0400694385A 6 - 0.0030161628A 7 + 0(A 8 ) (11) 

Table V shows estimates of the critical point A c and critical index from both the spin gap and susceptibility series 
for the coupling ratios J j J' = 1.25 and 1.40, using unbiased Dlog Pade approximants, where for J/J' = 1.40, we 
show the results for the spin gap at both k = and k = (7r/4, 37r/16). Apart from the spin gap A(7r/4, 37r/16) 
at J j J 1 = 1.40, there is considerable scatter in the results, but several features are apparent. The spin gap series 
show decreasing estimates of A c with increasing order for both J I J' ratios. The data are more consistent with the 
conclusion A c < 1, rather than A c > 1. For the spin gap A(7r/4,37r/16) at J j J' = 1.40, the unbiased Dlog Pade 
approximants give very convergent results, with the critical point A c = 0.988 and critical index v = 0.43. This 
indicates instability of the plaquette phase at A = 1, at least for this coupling ratio. This also indicates that there is 
indeed an intermediate phase. Note that the critical index v obtained here is about the same as that obtained from 



10 



the dimer series, suggesting that this transition might belong to the same new universality class. Our analysis shows 
that the transition probably does not belong in the classical 3D Hcisenberg universality class for J/ J' > 1.1. Because 
the critical index v for the minimum gap is much smaller than 1, we have obtained the dispersion results shown in 
Fig. u\ by performing the series extrapolation in a new variable 



s = i-(i-\/\ c y 



(12) 



in order to take the singular behavior into account. 

Finally we consider the second plaquette expansion (PE2). The minimum triplet gap is no longer at k = 0, at least 
for small to moderate A. This can be seen from the first few terms of the series 



A(k x ,k y )/J' = 1 - — (cosk x +cosk y ) + — cos(k x + k y ) 



(13) 



which gives the minimum gap at k x — k y = arccos( J' / J). It seems likely that the minimum triplet gap will remain at 
k^ 0, even for A ~ 1, although the series in this region are too erratic to confirm this. For this choice of unperturbed 
Hamiltonian we are also able to compute the excitation spectrum for the singlet excitation and we have obtained 
series to order A 7 . The minimum singlet gap is at k = and appears to be smaller than the triplet gap. An attempt 
to locate the critical point A c by Dlog Pade approximants to this series was hampered by poor convergence. We have 
attempted other analysis procedures, but without great success. But since the ground state energy from the second 
plaquette expansion (PE2) is higher than that obtained from both the first plaquette expansion and the Neel ordered 
Ising expansion, the probability of having the second plaquette configuration as the intermediate phase is remote. 




(b) J/J'=2 









FIG. 8. The triplet gap A(k)/J' versus A for k = (0, 0), (tt, 0) and (0, vr) and J/ J' = 0.875, and 2. Several different integrated 
differential approximants to the series are shown. 
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D. The energy gap in the columnar dimer expansion 

Finally we discuss the triplet dispersion obtained from the columnar dimer expansion. Here for J/ J' < 1 .15, we 
find the minimum gap is located at momentum k = (0,7r), as expected, since we expect to have a transition to the 
Neel ordered phase for small J/ J'. Fig. || show the gap at k = (0,0), (it, 0) and (0,7r) for J/J' = 0.875 versus A 
obtained from the integrated differential approximants |l2[ to the series. We can see that the minimum gap vanishes 
at a critical value A c < 1. We also expect this transition (to Neel order) to lie in the same universality class as the 
classical 3D Heisenberg model. To determine the phase boundary of columnar dimer order, we can use Dlog Pade 
approximants to both the series for the minimum gap and for the antiferromagnetic susceptibility. To get more reliable 
estimates for the critical point, we assume the critical exponents to be 7 ~ 1.4, v ~ 0.71. The results are shown in 
the phase diagram Fig. 0, where we can see that A c = 0.514(6) for J = 0, and that A c increases for increasing J . For 
J J J' — 1.15, we estimate A c = 0.97(5), but unfortunately the biased Dlog Pade results here are not accurate enough 
to tell whether A c is larger than 1. 

For J J J' > 1.75, we find that the minimum gap is located at k = (0, 0). Fig. g also shows the gap at k = (0, 0), 
(ir, 0) and (0,7r) for J/ J' = 2 versus A. Once again, one can locate the phase boundary by using the Dlog Pade 
approximants to the series for the minimum gap, and the results are also shown in Fig. |9j. We see that as J/ J' 
decreases, A c increases, but again, the analysis is not accurate enough to tell whether A c is larger than 1 in the 
intermediate region. 

In the analysis, we also note that the gap at k = (n, 0) shows some peculiar features, as shown in Fig. [|. For all 
values of J/J', the majority of integrated differential approximants to the series show that the gap increases for small 
A, then decreases dramatically to zero for A ~ 0.5. The majority of Dlog Pade approximants to the series also show 
that the series has a critical point at A ~ 0.5 with very small critical index (~ 0.2). But we believe this peculiar 
behavious is an artefact of the short series. 

For the most interesting region 1.15 < J/ J' < 1.75, the situation is quite complicated: the location of the minimum 
gap depends on both J/ J' and A. We take the mid-point for the presumed intermediate phase, J/ J' = 4/3, as 
example. For A near 1, we find the minimum gap for J/ J' — 4/3 is at k = (0.47, 2.80), slightly away from (0, 7r). The 
standard Dlog Pade approximants to the series for the minimum gap gives a critical point A c = 1.0(1), with critical 
index v ~ 0.4. This seems to be consistent with the indices obtained from the dimer and plaquette expansions, Since 
the standard Dlog Pade approximants cannot tell whether A c is larger than 1, we examine the extrapolations for the 
gap itself obtained from integrated differential approximants. The results for the gap at k = (0.47,2.80), (0,0) and 
(0, 7r) are shown in Fig. [10| The gap at A = 1 for k = (0.47, 2.80) is very small (~ 0.1J') when estimated from the 
direct A-series. When analysed in term of variable S (Eq. |l2|) the gap is compatible with zero. If the gap vanishes at 
or before A = 1, the columnar dimer phase would also be ruled out as intermediate phase for the Shastry-Suthcrland 
model. 
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FIG. 9. The phase diagram for columnar dimer order, where the curve in the small J/(J + J') region is determined by the 
Dlog Pade approximants to the gap at k = (0, n) and the antiferromagnetic susceptibility xaf, and we assume the transition 
lies in at the same universality class as the classical 3D Heisenberg model. The curve in the large J/(J + J') region is determined 
by the Dlog Pade approximants to the gap at k = (0, 0). 

IV. CONCLUSIONS 



We have attempted to further elucidate the nature of the phase diagram of the Shastry-Sutherland spin model at 
T = 0, by series expansion methods. We have significantly extended previous series, computed by ourselves J5| and 
others [|5|, and we have also derived a number of new series. The analysis of the various series allows us to draw some 
fairly firm conclusions, as well as others which are more tentative. 

The first question is whether an intermediate phase between the Nccl-ordercd phase and the dimer phase exists at 
all. It seems rather clear that the singlet dimer phase, with a simple singlet-product ground state, persists from large 
J I J' down to J J J 1 ~ 1.5 — 1.6. In our previous work |5J we argued that there is a direct first-order transition from 
the dimer phase to Neel order at this point. Our present results show the staggered magnetization and perpendicular 
susceptibility in the Neel phase vanishing at J/ J' ~ 1.2, with large uncertainty ±0.1. This appears more consistent 
with the Neel phase terminating at a second-order phase transition, so that there may indeed be an intermediate 
phase stretching over a very small range of coupling constants 1.20 < J/ J' < 1.5, in agreement with the suggestion 
of Koga and Kawakami B . 



Carpentier and Balents |16| have recently discussed an effective mean- field theory approach to the generalized 
Shastry-Sutherland model. They argue that a direct continuous transition from the Neel phase to the dimer phase 
is not possible. If this argument is accepted, it means that the vanishing of the order parameter and susceptibility 
in the Neel phase cannot be associated with a second order transition directly to the dimer phase. It follows that an 
intermediate phase must exist, over at least a small range of couplings. 

The next question is the nature of the transitions from the Neel and dimer phases, respectively, into the intermediate 
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FIG. 10. The triplet gap A(k)/J' versus A for k = (0, 0), (0.47, 2.80) and (0, tt) and J/ J' = 4/3. Several different integrated 
differential approxirnants to the series are shown. 



phase. The vanishing staggered magnetization and perpendicular susceptibility in the Neel phase indicate a 2nd order 
transition; but we are unable to estimate the critical exponents, since we do not have explicit series coefficients for 
an expansion in J/ J'. The vanishing energy gap would indicate that the transition from the dimer phase is also 
2nd order; but the result of Knetter et al. B showing that the two-particle bound-state energy vanishes before the 
single triplet gap indicates that more work needs to be done to characterize this transition also. It is interesting that 
one-particle gaps appears to vanish with a similar exponent in the dimer expansion, the plaquette expansion, and also 
the columnar dimer expansion. 

A further question concerns the nature of the intermediate phase. None of the suggested phases exhibit a ground- 
state energy which is distinctly lower than the Neel energy in the relevant coupling regime. Our results do not support 
the suggestion J3| that the intermediate phase is plaquette ordered. The extrapolated ground-state energy from the 
Neel expansion consistently lies below that from the plaquette expansion; and the triplet gap from the plaquette 
expansion appears to vanish at or before the physical value A = 1 is reached, indicating instability in this phase. Our 
longer series thus appear to contradict the conclusions of Koga and Kawakami |6| . Equally, our results do not support 
the suggestion that the intermediate phase is helically ordered [HM- The extrapolated ground-state energy from the 
Neel expansion consistently lies below that of the helically-ordered state, for any value of the spin orientation q other 
than 7T. We note that the mean-field approach of Albrecht and Mila 0] is strictly only valid for d > 2; while the 1/JV 
expansion approach of Chung, Marston, and Sachdev |7fl is of doubtful validity for the spin-j case, and is certainly 
not quantitatively accurate. 

We have also explored the possibility that the intermediate phase is a columnar dimer phase, as in the J\ — Ji 



Heisenberg model 14 , or the spin-i Heisenberg model on an anisotropic triangular lattice pa]. The extrapolated 
ground state energy for this state is comparable with that of the Neel state, within error bars, and the triplet gap from 
this expansion shows a window 1.15 ^ J/ J' < 1.5 in which it may remain finite, though very small, at the physical 



14 



value A = 1. This remains a possible candidate for the intermediate phase, although our data are not sufficiently 
accurate to allow a definitive conclusion. 

The nature of an intermediate state in this model thus remains an open question. It could, for example, be a 
structureless spin liquid, in which case a signature of this phase is elusive. Carpentier and Balents p6[ have also 
discussed possible intermediate phases in the generalized model, including a WISDQ (weakly incommensurate spin- 
density wave) ordered state -i.e. a periodic modulation of the expectation value of the total spin on each dimer - and 
a 'fractionalized' state, with topological order and deconfined spin-j excitations ('spinons'). We currently have no 
information to offer regarding these states. The fact that the energy gap drops to zero before or near the physical 
value A = 1 in every expansion we have tried appears to suggest that the intermediate phase is gapless, or near to 
it. Further support for this idea is given by the strong finite-size dependence of the exact diagonalization results of 
Miyahara and Ueda ||] in this region. 

The material SrCu2 (1303)2 appears to lie within the dimer phase of the Shastry-Sutherland model. Original 
estimates J/J' = 1.47 |^jj| P u t the material close to the phase transition point, and it was argued that some of the 
unusual properties of the material were due to this closeness. A more recent estimate J8|], including interactions in the 
third direction, gives J/ J' ~ 1.66, somewhat beyond even the point at which the two-particle gap vanishes. This puts 
the material clearly within the dimer phase and hence the question of intermediate phases is primarily of theoretical 
interest. 

It is clear that much work remains to be done on this fascinating model. It would be interesting to confirm and 
extend the results of Knetter et al. B concerning the transition from the dimer phase; and also to try and construct 
series coefficients for an expansion in J / J' in the Neel phase. New numerical methods are badly needed to explore 
the intermediate phase also: the best nmethod to employ here remains a puzzle. 
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TABLE I. Series coefficients for the ground-state energy per site Eq/(NJ') obtained from the Ising expansion about 
Neel-order (Neel), the helical-order (helical) expansion, the columnar dimer expansion, and from the plaquette expansion 
without diagonal bonds (PEf), and with diagonal bonds (PE2) for J/ J' = f .25, 1.4. Series coefficients of A 1 are listed. 



PEl 



PE2 



Helical 



Noel 



columnar dimer 



-5.000000000000 x 10" 1 
0.000000000000 
-3.575303819444xl0" 2 
-9.498878761574xl0" 3 
-6. 520377148279 XlO" 3 
-3.110192567684xl0" 3 
-2.051453756690X10" 3 
-1.134481706160X10" 3 
-7.196107350240X10" 4 



-4.218750000000X10" 1 

0.000000000000 
-1.577687449269X10" 1 

2.403841825921 xlO" 2 
-2.785619169517xl0" 3 

1.333554456619X10" 3 
-9.330263748650X10" 3 

9.793071057074xl0" 3 



J/J' = 1.25 
-3.562500000000x10 

0.000000000000 
-2.054711318598x10 

4.163537328971x10 
-4.669142343282x10 

2.832044653274x10 
-3.363486912444x10 

3.693336765114x10 
-5.318671447970x10 

7.567434418701x10 
-1.138657557136x10 

1.726680824433x10 



-3.437500000000x10 

0.000000000000 
-2.857142857143x10 

1.020408163265x10 
-9.273368951135x10 

1.575551772772x10 
-2.473894356829x10 

4.097681140990x10 
-8.162334082025x10" 

1.693646904376 
-3.494392524807 

7.456679516027 
-1.646408487950X10 1 



-3.750000000000x10 
0.000000000000 

-1.904296875000x10 
7.141113281250x10 
9.905815124512x10 

-9.499887625376x10 

-1.780716847214x10 
1.494069254463x10 

-3.803820899197x10 
9.729690462284x10 

-2.666137442909x10 



-5.000000000000 x 10" 1 
0.000000000000 
-3.392361111111X10" 2 
-8. 058449074074 xlO" 3 
-5.431377864048 xl0~ 3 
-2.536510876037xl0" 3 
-1.754982806912 xl0~ 3 
-1.017310795905X10" 3 
-7.033788209007x10" 



1-4 



-4.125000000000x10" 

0.000000000000 
-1.678565705128x10" 

2.952951931067x10" 
-1.145691038736x10" 

2.724731917375x10" 
-1.315792108890x10" 

1.282279331689x10" 



J/J' = 1.40 
-3.535714285714X10" 1 

0.000000000000 
-1.943166023166 xlO" 1 

3.629865386622X10" 2 
-4.218721381575X10" 2 

2.276170842636 xlO" 2 
-3.202093475641 x 10" 2 

2.995390698796 xlO" 2 
-4.311390066885x10 

5.834414925135x10 
-8.797286992572x10 

1.341535224738x10 



-3.250000000000x10" 

0.000000000000 
-3.125000000000x10" 

1.367187500000x10" 
-1.664227322346x10" 

3.302972722063x10" 
-6.320111458811x10" 

1.313875979565 
-3.190631154651 

8.019181935184 
-2.044153224739 xlO 1 

5.383666306276 xlO 1 
-1.457331428464X10 2 



-3.750000000000 x 10" 1 

0.000000000000 
-1.950000000000X10" 1 

1.659375000000X10" 2 

1.704765625000X10" 2 
-1.338582682292X10" 2 
-1.870258387587X10" 2 

1.669764147147X10" 2 
-1.134220766188X10" 3 
-6.547755059723X10" 3 
-1.637369775362 xlO" 3 
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TABLE II. Series coefficients for 
expansion without diagonal bonds 
listed. 



the triplet gap A/ J' at k = and staggered susceptibility xaf 
(PE1), and with diagonal bonds (PE2) for J/ J' = 1.25, 1.4. 



obtained from the plaquette 
Series coefficients of A 1 are 



PEl 



PE2 



J/J' = 1.25 



J I J' = 1.4 



J I J' = 1.25 



J I J' = 1.4 







triplet gap A/./' 


at k = 









1.000000000000 


1.000000000000 




1.000000000000 


1.000000000000 


1 


-5.000000000000X10 -1 


-4.000000000000X10 -1 




-9.166666666667xl0~ 1 


-8.666666666667x10" 


2 


-1.949508101852xl0~ 1 


-1.558796296296xl0~ 1 




5.848101551222xl0~ 4 


-1.753739316239x10" 


3 


-1.215547789271 xl0~ 2 


6.279063786008xl0~ 3 




3. 919269654832 xlO -1 


7.809106177672x10" 


4 


-4.262073088443 xl0~ 2 


-3.063595079003xl0~ 2 




-1.209759818705 


-2.966617086045 


5 


-1.169575496893xl0~ 2 


-4.837854449338X 10~ 3 




2.808227188528 


1.165893864520X10 1 


6 


-2.087090156347xl0~ 2 


-1.752776610250xl0~ 2 




-7.118402535782 


-5.396727727606 xlO 1 


7 


-9.640627326157xl0~ 3 


-8.067626260181 xl0~ 3 









1.333333333333 
1.333333333333 
1.250353652263 
1.201753927683 
1.113460517630 
1.044513939707 
9.707650534281x10" 



staggered susceptibility xaf 
1.333333333333 
1.066666666667 
7.572093621399 xlO -1 
5.539876821845X10 -1 
3. 809246267771 x 10 _1 
2.747370398592 xlO -1 
1.964730571491xl0~ 1 



1.333333333333 
2.444444444444 
2.323016302459 
2.280972931735 
2.075318140109 
2.415891359693 
2.233378054361 



1.333333333333 

2.311111111111 

1.684709846149 

1.010069843149 

3.056543491121x10 

6.618315713507x10 

4.157664670338x10 



TABLE III. Series coefficients for the stag, 
and the helical-order (helical) expansion for 



;gered magnetization M obtained from the Ising expansion about Neel-order (Neel) , 
J I J 1 — 1.25, 1.4. Series coefficients of \ l are listed. 







helical 






Neel 


i 


J/ J' = 1.25 




J/ J' = 1.4 




J/ J' = 1.25 


J/ J' = 1.4 





5.000000000000x10" 


-l 


5.000000000000x10" 


-l 


5. 000000000000 xlO -1 


5.000000000000x10" 


1 


0.000000000000 




0.000000000000 




0.000000000000 


0.000000000000 


2 


-1.976862628853x10" 


-l 


-1.787535427319x10" 


-l 


-3.265306122449 x 10" 1 


-3.906250000000x10" 


3 


8.015745892489x10" 


-2 


6.696385209168x10" 


-2 


2.332361516035xl0~ 1 


3.417968750000x10" 


4 


-2.020938451061x10" 


-1 


-1.703648728139x10" 


-1 


-6. 249672540570 xlO -1 


-1.188130158850 


5 


1.595089057636x10" 


-1 


1.186622673353x10" 


-1 


1.247039703760 


2.891792139388 


6 


-2.856270160505x10" 


-1 


-2.387988493110x10" 


-1 


-2.564444941065 


-7.392698990572 


7 


3.509646624490x10" 


-1 


2.550333854825x10" 


-1 


5.414869562182 


1.960927054837X10 1 


8 


-6.159923626864x10" 


-1 


-4.541054989763x10" 


-1 


-1.300042701432 xlO 1 


-5.777973112981 xlO 1 


9 


9.506010225594x10" 


-1 


6.635675648444x10" 


-1 


3.100311343644X10 1 


1.688382921512X10 2 


10 


-1.636489163183 




-1.143031741464 




-7.335660790958 XlO 1 


-4. 957990614738 x 10 2 


11 


2.726180882130 




1.887668026020 




1.769340004349 XlO 2 


1.478774628135 xlO 3 


12 










-4.347120802018X10 2 


-4.471290404091 xlO 3 



TABLE IV. The results of {m/n/l} 
J I J' = 1.4 from the Neel ordered Ising 
denotes a defective approximant. 



integrated differential approximants to the series for the ground state Eo/NJ' for 
expansion and the plaquette expansion without diagonal bonds (PE1). An asterisk 



n 


{{n-3)/n/l} 


{{n-2)/n/l} 


{(n - l)/n/l} 


{n/n/l} 


{{n + \)/n/l} 


{{n + 2)/n/l} 


{{n + 3)/n/l} 








Nccl ordered 


Ising expansion 
= 

-0.54979 








n= 3 




-0.52170 


-0.55795 


-0.57974 


-0.55498 


-0.57207 


n= 4 


* 


-0.55178 


-0.55573 


-0.56195 


-0.56468 


-0.56026 


-0.55991 


n= 5 


-0.56212 


-0.57002 


-0.56907 


-0.56264 


-0.55989 






n= 6 


-0.56912 


* 


-0.56121 










n= 7 


-0.56405 




/ 


= 1 

-0.55765 








n= 3 




-0.55300 


-0.55127 


-0.57441 


-0.56242 


-0.55104 


n= 4 


-0.55139 


* 


-0.55545 


-0.56380 


* 


-0.56002 




n= 5 


-0.56139 


-0.56960 


-0.56066 


-0.56100 








n= 6 


-0.56503 


-0.56097 


/ 


= 2 

-0.55148 








n= 2 




-0.56352 


* 


-0.55131 


-0.60526 


-0.57327 


n= 3 


-0.55123 


-0.55086 


-0.55252 


* 


-0.56628 


-0.54598 


-0.55427 


n= 4 


-0.55270 


* 


-0.55923 


-0.55915 


* 






n= 5 


-0.56600 


-0.55917 


-0.56079 










n= 6 


-0.56213 




/ 


= 3 

-0.55914 








n= 1 






i 
-0.56410 


-0.55010 


-0.56012 


-0.57599 


n= 2 




-0.56089 


* 


-0.55363 


* 


-0.56992 


* 


n= 3 


-0.55104 


-0.55440 


-0.55944 


-0.55918 


-0.56186 


-0.56198 




n= 4 


* 


-0.55918 


-0.55923 


-0.56196 








n= 5 


-0.56143 


-0.56258 













plaquette expansion without diagonal bond (PE1) 
I = 



n= 1 










-0.55534 


-0.55281 


-0.55647 


n= 2 






-0.55246 


-0.55301 


-0.55390 


-0.55473 


-0.55482 


n= 3 




-0.55297 


* 


-0.55492 


-0.55484 






n= 4 


-0.55425 


-0.55448 


-0.55484 










n= 5 


* 






I = 1 








n= 1 








* 


-0.55243 


* 


-0.55344 


n= 2 






-0.55316 


-0.55574 


-0.55595 


-0.55499 




n= 3 




* 


-0.55587 


-0.55467 








n= 4 


-0.55458 


-0.55513 




I = 2 








n= 1 






* 


-0.55270 


* 


* 


* 


n= 2 




-0.55347 


-0.55481 


* 


-0.55488 






n= 3 


* 


-0.55496 


-0.55487 










n= 4 


-0.55492 
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TABLE V. [n/m] Dlog Pade approximants to the series for the triplet gap A(k) and antifcrromagnetic susceptibility xaf 
for the plaquette expansion without diagonal bonds (PE1). An asterisk denotes a defective approximant. 



n 


[(« - 2)/n] 


[(« - l)/«] 


[n/n] 


[(«+!)/«] 


[(n + 2)/n] 




pole (residue) 


pole (residue) 


pole (residue) 


pole (residue) 


pole (residue) 








A(0, 0) for J/ J' = 1.25 






n= 1 






1.4098(1.2718) 


0.8592(0.2879) 


1.2001(1.0957) 


n= 2 




1.0804(0.6523) 


1.0530(0.6018) 


1.0305(0.5502) 


1.0181(0.5170) 


n= 3 


1.0555(0.6075) 


0.9247(0.2175)* 


1.0074(0.4801) 






n= 4 


1.0089(0.4855) 
















A(0,0) for J/J' = 1.4 






n= 1 






2.0315(1.9470) 


0.8107(0.1237) 


1.5546(1.6732) 


n= 2 




1.2949(0.6512) 


1.2016(0.5172) 


1.1075(0.3712) 


1.0273(0.2532) 


n= 3 


1.2181(0.5462) 


0.7442(0.0227)* 


0.9352(0.1267) 






n= 4 


0.9454(0.1386) 














A| 


>/4,37r/16) for J/ J' = 1.4 






n= 1 






1.6055(1.4266) 


0.638(0.0896) 


1.4385(2.3132) 


n= 2 


* 


0.9883(0.4300) 


0.9889(0.4307) 


0.985(0.4240) 


1.0203(0.5055) 


n= 3 


0.9889(0.4307) 


0.9884(0.4301) 


0.9885(0.4302) 






n= 4 


0.9885(0.4302) 
















XAF for J/ J' = 1.25 






n= 1 






0.9830(-0.8461) 


1.2245(-1.6351) 


0.9962(-0.7165) 


n= 2 


1.0003(-0.8898) 


1.0739(-1.0477) 


1.1011(-1.1273) 


1.0825(-1.0510) 




n= 3 


1.1033(-1.1367) 


1.0910(-1.0902) 
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